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SUMMARY 

An  approximate method  is  presented  for  calculating  the  pressure 
distribution  on  conical  bodies  of  noncircular  cross  section  in  a  super¬ 
sonic  flow  field.  By  a  superposition  of  elementary  conical  flows  due 
to  line  sources,  the  flow  about  an  arbitrary  cone  may  be  described.^ 
Illustrations  of  the  pressure  distribution  about  several  shapes  are 
included  to  demonstrate  the  method.  The  problem  of  such  a  body  at 
angle  of  attack  may  also  be  solved  by  the  same  method,  as  well  as  the 
problem  of  yawed  flight  at  angle  of  attack. 


INTRODUCTION 

Linearized  methods  have  been  applied  to  determine  solutions  of 
the  supersonic  flow  field  about  both  solid  and  open-nosed  bodies  of 
revolution.  (See  references  1  'to  4.)  The  methods  of  reference  1  have 
been  extended  herein  to  produce  a  linearized  solution  for  cones  of  arbi¬ 
trary  cross  section,  such  as  might  serve  as  forebodies  of  nonsymmetrical 
fuselages.  The  solution  is  based  on  the  use  of  a  combination  of  line 
sources  inclined  arbitrarily  to  the  flow  direction.  If  the  proper  sources 
are  chosen,  any  conical  body  shape  can  be  described  and  the  resulting 
surface  pressures  can  be  calculated. 

The  method  presented  herein  was  devised  during  the  fall  of  1947 
at  the  KACA  Cleveland  laboratory. 


SYMBOLS 

The  following  symbols  are  used  in  this  analysis: 

C  pressure  coefficient 
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constant  proportional  to  source  strength 
free -stream  Mach  number 

slope  of  line  source  with  respect  to  flow  direction 

component,  normal  to  flow  direction,  of  distance  from  point 
'on  line  source  to  arbitrary  point  in  flow  field 

free-stream  velocity 

component  of  total  velocity  normal  to  body  at  surface 

radial  (cylindrical  coordinate)  perturbation-velocity  component 

axial  perturbation-velocity  component 

tangential  (cylindrical  coordinate)  perturbation  velocity  component 
cylindrical  coordinates 
cotangent  of  Mach  angle,  V  - 1 

angular  position  of  line  source,  measiired  from  6  =  tt/2  plane 

angle  between  radial  velocity ^and  normal  to  body,  measured  in 
X  =  constant  plane 

angle  whose  tangent  is  ratio  of  radial  to  axial  coordinates 
of  point  on  body  surface,  tan”^  r/x 

axial  coordinate  of  line  source 

source  strength  per  unit  of  axial  length 

perturbation-velocity  potential 

angle  between  normal  to  body  surface  and  normal  to  ray, 
which  intersects  axis  of  body 


GEIJERAL  ANALYSIS 

In  the  solution  of  the  pressure  distribution  over  a  body,  a  method 
of  successive  approximation  is  used*  A  perturbation-velocity  potential 
based  on  linearized  flow  is  found  from  which  the  three  velocity  components 
may  be  determined-  If  these  velocities  satisfy  the  boundary  conditions 
for  the  desired  body,  the  potential  describes  the  flow  about  the  actual 
body.  The  pressure  coefficient  may  then  be  found  from 
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C  =  -  2U  /U 
P  X 


(1) 


potential  may  not  be  assumed  completely 
tut^is  subject  to  several  general  limitations  and  to  some 
^rticular  ones  imposed  by  the  body.  The  particular  limitations  form 

arf  potential.  Ihe  two  general  limitations 

JIV  potential  must  satisfy  the.Prandtl-Glauert  equation;  and 

of  linearized  flow,  the  disturbanL  velocities 
should  disappear  at  the  Mach  cone. 

■nrvT.  Prandtl-G^uert  equation  for  compressible  frictionless  potential 
flow^  in  linearized  form  and  cylindrical  coordinates  is 


(m2  .  1) 


bx  r  dr 


dr 


O  CD 

o  - 5  =  0 

r^  de 


(2) 


•  source  or  of  strength  per  unit  length  f(nj  lying 

in  the  plane  (.  -  [(tt/s)  +  6  J,  whose  slope  with  respect  to  the ^f low 

direction  is  m,  is  shown  in  figure  1.  The  potential  at  seme  point  P 
due  to  such  a  disturbance  is  i'oinu  r 


X-R{3 


f  (1)  d  e 


(3) 


'0  aJ  (  I  - 

where  R  Is^  from  the  geometry  of  figure  1,  defined  by 

r2  =  r^  +  (jii|)2  _  2infr  sin  (6  -  6)  (4) 

^mtloJ^(2r^  readily  be  shown  by  direct  substitution  to  satisfy 


in 


Combining  equations  (3)  and  (4)  and  rearranging  the  terms  res;at 
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’rom  figure  2(a), 

U  =  0  =  U  cos  i?  -  (U  +  U  )  sin  V  cos  sin  \i/ 

nr  X  y 


Ur  -  Ug 


— =  (U  +  Uy)  tan  V 

COS  V  '  x^ 


But,  also  from  figure  2(a), 


tan  ^  =  tan  X  cos  V 


and,  with  reference  to  figure  2(1), 


tan  X  = 


where  r  =  f(0)  is  the  equation  of  the  hody  cross  section.  If  equa¬ 
tions  (lO)  to  (12)  are  combined,  the  boundary  condition  for  the  flow 
is  obtained: 


u  -  i  ^  U 
r  de 


-  ('j  +  ''*)  I 


The  second  general  condition  for  linearized  conical  flow,  stated 
previously,  is  that  the  perturbation  velocities  approach  zero  as  the 
Mach  cone  is  approached,  or  as  rP/x  approaches  imity.  This  con¬ 
dition  is  exactly  true  for  equations  (7)  to  (9). 

The  potential  found  from  the  single-line  source  that  has  been 
considered  until  now  is  Insufficient  to  calculate  the  flow  about  an 
arbitrary  body.  Because  the  potential  due  to  one  source  satisfies 
the  general  limitations  of  the  problem,  however,  a  series  of  poten¬ 
tials  due  to  a  number  of  sources  of  various  strengths  and  positions 
can  be  so  added  together  that  the  resulting  flow  satisfies  the  bound¬ 
ary  conditions  for  the  particular  body  in  question  and  thus  can  be 
assumed  to  be  the  desired  flow. 

Foiir  parameters  are  considered  in  selecting  the  source  pattern: 
the  angular  position  6,  the  slope  of  each  source  relative  to  the 
flow  direction  m,  the  number  of  sources,  and  the  strength  of  each. 


APPLICATION  OF  METHOD 


Because  the  perturbation  velocities  are  fimctions  of  only  r/x 
and  6,  the  resulting  body  contours  are  conical  and  only  one  body 


section  that  is  in  a  plane  normal  to  the  free-streem-flov  direction 
■need  he  considered.  The  problem  is  thereby  reduced  to  two  dimensions 
and  the  body  is  completely  described  by  one  plane,  as  in  figure  3. 

The  line  sources  thus  appear  in  such  a  figure  simply  as  points. 

With  the  use  of  several  sources  of  varying  strength  and  position, 
the  boundary  condition,  defined  by  equation  (13) ,  can  be  made  to  yield 
the  flow  about  the  desired  body  shape.  Although  such  a  solution  is 
found  largely  by  trial  and  error,  some  general  rules  for  the  selection 
of  sources  yielding  a  specified  body  sliape  can  be  established.  The 
four  parameters  previously  mentioned,  that  is,  the  angular  posi¬ 
tion  5,  the  slope  of  the  line  sources  m,  the  number  of  sources,  and 
the  source  strength,  must  be  kept  in  mind. 

1.  The  axes  of  symmetry  of  the  body  section  should  be  noted, 
inasimach  as  the  sources  must  be  symmetrically  arranged  relative  to 
the  same  axes. 

2.  The  source  nearest  to  a  peak  in  the  section  should  be  nearer 
to  tliat  point  than  to  any  other  point  on  the  section  because  the  per¬ 
turbation  velocities  due  to  the  body  are  a  maximum  at  the  peak  and, 
inasmuch  as  velocities  due  to  the  source  increase  as  the  distance 
from  the  source  decreases,  the  point  closest  to  the  source  has  the 
hipest  velocity.  The  distance  from  a  point  on  the  body  to  the 
nearest  source  must  therefore  be  equal  to  or  less  than  the  radius 

of  curvature  of  the  body  section  at  that  point. 

5.  If  the  body  is  elongated,  a  series  of  sovarces  in  a  line  are 
required.  The  slenderer  the  section,  the  more  sources  are  needed 
to  prevent  contour  irregularities  in  the  body  described  by  the 
resulting  solution  (equation  (15)).  Also,  for  a  given  body  of  this 
type,  as  the  section  najrrows  the  sources  should  be  closer  together. 

4.  In  general,  the  larger  the  number  of  sources  used,  the  more 
closely  the  linearized  flow* obtained  fits  the  body  in  question  and 
the  fewer  are  the  number  of  trial  solutions  required  to  obtain  a 
satisfactory  answer.  Each  such  solution,  however,  is  more  laborious 
than  one  using  fewer  sources. 

To  summarize,  the  angular  position  of  the  sources  is  determined 
by  rules  1  and  2,  the  distance  from  the ^axis  by  1,  2,  and  3,  and  the 
number  of  sources  by  1  and  3.  Rule  4  serves  as  an.  over-all  guide. 

Only  K,  which  is  proportional  to  the  source  strengths,  remains  to 
be  determined.  Erom  the  symmetry  condition,  sources  tliat  are  in  the 
same  position  relative  to  the  axes  of  symmetry  have  the  same  strength. 
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With  a  known  number  of  different  source  strengths^  equation  (13)  may 
be  solved  for  the  strengths  at  the  same  number  of  points  on  the  surface 
by  using  equations  (?)  to  (9).  The  potential  is  now  completely  defined. 
This  potential  shoxild  then  be  checked  in  equation  (13)  at  several  points 
on  the  body  to  determine  whether  the  flow  due  to  the  soixrce  potential 
is  the  same  as  that  over  the  actual  body.  It  should  be  noted  that 
the  strength  of  a  source  may  be  negative. 

The  case  of  flow  at  angle  of  attack  can  be  solved  by  considering, 
at  zero  angle,  a  body  whose  cross  section  consists  of  sections  of  the 
actual  body  taken  normal  to  the  free  stream  instead  of  normal  to  the 
axis.  The  x-axis  of  the  new  body  is  then  parallel  to  the  flow  direc¬ 
tion.  For  a  small  angle  of  attack,  the  two  cross  sections  probably 
differ  little,  but  their  positions  relative  to  the  coordinate  axes 
differ.  The  case  of  yawed  flight  at  angle  of  attack  can  be  solved 
in  a  similar  manner  (fig.  4). 


Examples 

Several  examples  follow  to  illustrate  the  general  rules  that 
have  been  outlined: 

Example  I.  -  Assume  that  an  elliptic  section  (fig.  3).  is  desired. 
From  the  symmetry  condition,  two  sources  are  assumed  as  a  first  approx¬ 
imation  to  such  a  section.  These  sources  are  of  equal  strength,  equi¬ 
distant  from  the  x-axis,  and  at  5=0  and  6  =  n.  The  velocities 
are  found  from  equations  (7)  to  (9). 

The  equation  of  an  ellipse  is 

2 

p  =  — - - - -  - - -  -  . . I  ■  .  I 

2  2  2  2 
d  sin  0  +  c  cos  B 

where  c  and  d  are  major  and  minor  semi -axes,  respectively. 

1  dr 

From  this  relation,  -  ^  may  be  obtained: 


1  ^  _  sin  B  cos  0'(c^-d^) 

^  sln^  e  +  cos^  e 


(14) 


\Ttien  a  Mach  number  and  values  for  c  and  d,  are  given,  a  position 
for  the  sources  (that  is,  a  value  of  m)  is  selected.  From  the 
second  general  rule,  c  -  m  of  figure  3  should  be  somewhat  less  than 


•the  radius  of  curvature  at  c,  n/2.  With  this  guide,  a  value  of  m 
can  he  assumed..  From  the  coordinates  of  one  point  on  the  surface  K/U 
can  he  found  hy  substitution  of  equations  (7)  to  (9)  and  (14)  in  equa¬ 
tion  (13)*  Several  other  points  should  then  he  checked  in  the' boundary 
relation.  If  agreement  is  poor,  the  sources  should  he  moved  and  per- 
iiaps  others  added.  When  a  potential  is  obtained  that  satisfies  equa¬ 
tion  (13)  closely  enough,  the  pressure  coefficient  can  then  he  found 
from  equation  (l). 

An  example  of  the  results  of  such  a  calculation  Is  shown  in 
figure  3.  The  body  shape  obtained  by  using  the  value  of  K  previ¬ 
ously  determined  and  by  solving  equation  (13)  for  r/x  as  a  function 
of  6  is  shown  together  with  the  desired  shape.  In  this  example,  where 
the  contour  is  almost  circular,  the  pressure  distribution  approximately 
follows  the  body  shape.  The  deviation  of  the  calculated  section  from 
the  desired  ellipse  could  be  considerably  decreased  by  placing  the  two 
sources  farther  apart  and  adding  a  third  source  at  the  origin. 

If  the  flow  about  an  ellipse  having  a  larger  ratio  of  major  to 
minor  axes  than  was  used  in  figure  3  is  desired,  the  perturbation-velocity 
potential  for  two  sources  gives  poor  results.  The  effect  of  varying'  K 
while  holding  m  constant  is  shown  in  figure  5.  A  similar  result  is 
obtained  by  holding  the  length  of  one  of  the  axes  constant  and  varying 
the  length  of  the  other  by  changing  m,  while  holding  the  source  strength 
constant.  In  order  to  obtain  a  satisfactory  solution  for  an  ellipse  having 
the  ratio  of  the  axes  much  greater  than  that  in  figure  3,  a  series  of  sources 
can  be  used. 

Example  II.  -  Assiune  that  an  ellipse  having  the  ratio  of  major  to 
minor  axes  equal  to  3  is  desired.  (See  fig.  6.)  Obviously,  for  symmetry, 
a  number  of  source  pairs,  as  used  in  the  previous  example,  plus  perhaps 
a  single  source  at  the  origin,  will  give  the  desired  solution.  The  position 
of  source  1  in  figure  6  may  he  assumed  from  the  condition  that  its  distance 
from  the  peak  must  be  approximately  the  radius  of  curvature  at  the  peak. 

The  distance  between  sources  1  and  2  is  taken  as  about  equal  to  the  distance 
from  source  2  to  the  nearest  point  on  the  body.  The  remaining  sources  may 
be  similarly  chosen.  This  procedure  gives  a  system  of  three  source  pairs 
whose  velocities  are  found  from  equations  (?)  to  (9),  and  a  single  source 
whose  velocities  are  found  from  equations  (7a)  and  (8a),  which  are  simply 
the  velocities  found  when  m  is  zero.  When  equation  (13)  is  solved  for 
this  system  of  sources  at  four  points  on  the  body,  K^,  E^,  K3,  and  E4 
(fig.  6^  are  determined  and  the  pressure  coefficient  can  be  found  from 
equation  (l). 

The  desired  ellipse  and  the  contour*  calculated  by  the  use  of  seven 
sources  are  shown  in  figure  6  together  with  the  pressure  distribution 


corresponding  to  the  source  configuration.  In  this  example  the  section 
deviates  greatly  from  a  circle  and  the  pressure  distribution  no  longer 
follows  the  body'  contour.  The  maximum  and  minimum  values  of  the  pres¬ 
sure  coefficient,  however,  still  occur  at  the  maximum  and  minimum  points, 
respectively,  on  the  section. 

Example  III.  -  Now  assiune  that  the  pressure  distribution  is  desired 
over  a  body  whose  cross  section  is  a  triangle  modified  by  rounding  the 
vertices  (fig.  7).  From  symmetry  considerations,  three  soinrces  of  equal 
strength  should  be  assumed  at  6=0,  2jt/3,  and  4jt/3.  These  sources 

should  be  placed  at  approximately  the  center  of  curvature  of  the  vertices. 
Because  the  body  is  approximately  circular,  a  source  that  is  not  of  the 
same  strength  as  the  others  should  be  placed  at  the  origin. 

The  velocity  components  can  then  be  calculated  from  eqimtlons  (7) 
to  (S)  and  the  strength  of  the  sources  can  be  established  by  solving 
equation  (13)  at  two  points. 

Such  a  surface,  with  the  corresponding  pressure  distribution.  Is 
Illustrated  in  figure  7.  As  in  the  first  example,  the  pressure  distri¬ 
bution  follows  the  trends  of  the.  body  shape. 

It  must  be  remembered  that  the  examples  given  are  meant  to  illustrate 
the  method  of  solution  of  such  bodies  rather  than  to  show  actual  pressure 
distributions,  although  the  trends  indicated  should  be  correct.  The  bodies 
chosen  are  probably  not  slender  enough  for  great  accuracy  in  a  linearized 
solution. 


SUMMAEY  OF  ANALYSIS 

An  approxlraatem  method  has  been  presented  for  calculating  the 
pressure  distribution  on  conical  bodies  of  arbitrary  cross  section  in 
supersonic  flow.  By  a  combination  of  elementary  conical  flows  due,  to 
line  sources,  the  flow  about  a  slender  arbitrary  cone  can  be  described. 
Four  parameters  are  considered  in  determining  such  a  system  of  sources : 
the  spacing  of  the  sources  around  an  axis  lying  in  the  flow  direction, 
the  slope  relative  to  the  flow,  the  niamber  of  sources,  and  their  strength. 
The  first  three  parameters  can  be  determined  by  several  conditions. 

First,  the  same  symmetries  will  hold  for  the  sources  as  for  the  body. 
Second,  the  distance  from  a  peak  to  a  source  will  be  less  than  or  equal 
to  the  radius  of  curvature  at  that  po'int.  ^ird,  for  an  elongated  body, 
the  slenderer  the  section,  the  closer  together  the  sources  must  be. 
Finally,  the  greater  the  number  of  sources  used,  the  more  accurately 
the  desired  body  can  be  approximated.  The  solution  will  be  more  labo¬ 
rious,  however,  if  more  sources  are  used.  The  fourth  parameter,  source 
strength,  may  be  found  by  direct  calculation. 
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§ 

The  case  of  a  conical  body  at  angle  of  attack  and  of  yawed  flight 
^  can  be  solved  by  the  same  method. 

Flight  Propulsion  Research  Laboratory > 

rfetional  Advisory  Committee  for  Aeronautic^, 

Cleveland,  Ohio,  April  27,  1948. 
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(a)  Projection  sketch  of  body  surface  and  velocity  components. 


(b)  Section  normal  to  x-axls  taken  through  point  0\ 

Figure  2.  -  Geometric  relations  between  body  surface  and  velocity 

components. 
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. .  Desired  section 

-  Calculated  section 

- Pressure  distribution 

O  Source  position 


m 


Figure  3,  -  Comparison  of  desired 
resulting  pressure  distribution 


elliptical  body  with  calculated  body  and 
.  p  =  1.5;  m  =  0.125;  K/U  =  0.0220, 
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Plgiire  4.-  Comparison  of  body  at  zero  angle  of  attack 

and  zero  yaw  with  same  body  at  angle  of  attack  and  In 
yawed  flight. 


Figure  5.  -  Surfaces  obtained  with  sources  of  varying  strength 

at  (a,  i  it/S).  m  =  0.20. 
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Kl/U  *  0.00483 
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*=  .19 

Kg/U  «  .00395 

- "“  —  Calculated  body 

®3 
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Figure  6,  -  Comparison  of  desired  slender  elliptical  body  with  calculated  body  and  resulting 

pressure  distribution.  p  =  1.5. 
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Desired  and  calculated  body 
Pressure  distribution 
O  Source  position 


Figure  7,  -  Comparison  of  desired  modified  triangular  body  with  calculated 
body  and  resulting  pressure  distribution.  {3  =  1,5;  m  =  0,16;  Kt/U  =  0,0102; 
Kg/U  =  0.0154. 
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